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ABSTRACT 

In  this  paper -we- study- the  system  (1.1),  (1.3),  which  describes  the 
stationary  motion  of  a  given  amount  of  a  compressible  heat  conducting,  viscous 

is*  ,..  ,  t  » 

fluid  in  a  bounded  domain  ft  of  rf1,  n  >  2.  Here  u(x)  is  the  velocity 
field,  p(x)  is  the  density  of  the  fluid,  £(x)  is  the  absolute  temperature, 

y 

f(x)  and  h(x)  are  the  assigned  external  force  field  and  heat  sources  per 
unit  mass,  and  p(p,£)  is  the  pressure.  In  the  physically  significant  case 
one  has  g  ■  0.  We  prove  that  for  small  data  (f ,g,h)  there  exists  a  unique 
solution  (u,p,C)  of f problem  (1.1),  (1.3)^,  in  a  neighborhood  of  (0,m,C0)> 
for  arbitrarily  large  data  the  stationary  solution  does  not  exist  in  general . 
(see  section  5) .  Moreover,  we  prove  that  (for  barotropic  flows)  the 
stationary  solution  of  the  Navier-Stokes  equations  (4.15)  is  the 
incompressible  limit  of  the  stationary  solutions  of  the  compressible  Navier- 

Stokes  equations  (4.14),  as  the  Mach  number  becomes  small.  Finally,  in 

. —  it  ^ 

section  5  we  will  study  the  equilibrium  solutions  for  system  (4.1).  For  a 

more  detailed  explanation  see  the  Introduction. 
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1.  Introduction 

In  this  paper  we  study  the  system 


(1.1) 


-uAu  -  v7  div  u  +  7p(p  ,C  )  ■  p  [f  -  (u*7  )u]  , 

div(pu)  -  g. 


-  X*C  +  cypu*7c  ♦  Cp£  (P  »C,)div  u  -  ph  +  i>(u,u),  in  0 


°'  5|r  -«o 


in  a  bounded  open  domain  0  in  eF1,  for  arbitrarily  large  n  >  2.  It  is  assumed  that 
0  lies  (locally)  on  one  side  of  its  boundary  T,  a  C2  manifold.  Here, 


*(u,u)  -  x 


n  3u.  3u,  2 

•  ,5^  *  ^ 


+  x ,(div  u) 


n 

and  (v*7 )u  -  J  v. (3u/3 x, ) .  System  (1.1)  describes  the  stationary  motion  of  a 
i“1 

compressible,  heat-conductive,  viscous  fluid  (see  Serrin  [10]).  In  equation  (1.1),  u(x)  - 
(u1(x),u2(x),...,un(x))  is  the  velocity  field,  p(x)  is  the  density  of  the  fluid,  C  (x) 
is  the  absolute  temperature,  f(x)  and  h(x)  are  the  assigned  external  force  field  and 
heat  sources  per  unit  mass,  and  p(P,C)  is  the  pressure.  In  the  physically  significant 
case  one  has  g  «  0,  however,  it  is  not  without  interest,  from  a  mathematical  point  of 
view,  to  study  the  general  case. 

In  order  to  avoid  technicalities,  we  will  assume  that  the  coefficients  v  >  0,  v  >  -y  , 
X  >  0,  cy,  Xg,  Xi»  are  constant.  A  dependence  of  those  coefficients  on  u,  P,  C,  as  well 
as  the  introduction  in  equation  (1.1)  of  other  hinds  of  non-linearities ,  does  not  give  rise 


to  substantial  difficulties.  For  the  same  reason,  we  assume  that  Cn  >  0  is  constant. 
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p(P,C)  £  C3+l(I(l,ij)).  _If_  (f,g,h)  e  wi+1'P  X  w^+2'p  x  w^'P  verify  the  condition 


Ifl.^.  +  ♦  Ihl  <  o’ 

J+1#P  j+2,p  j+1,p  o 


then  there  exlata  a  unique  solution  (u,p,;)  e  W^+3,p  x  w^+2,p  x  w3+3,p  of  problem  (1.1), 
(1.3)  in  the  ball 


lul^,  *  Ip -ml  +  lc-5  I  <  c* 

j+3,p  j+2,p  s^o  j+3,p  1 


In  the  above  statement  c'0,  c ' 1  are  suitable  positive  constants  depending  only  on  0, 


n,  p,  j,  P,  v,  k,  m,  C0>  e^,  x*  X0»  1/  tj.  Ti#  S^d-1,2),  where  -  sup  |w  i  «J /*  )  |  for 

,2+1, 


(o,a)  «  1(1, 1^),  and  is  the  norm  of  u)^(o, a)  in  the  space  C  J(X(t,t^)). 

The  LP-theory  enables  us  to  treat  the  n-diaensional  case  by  handling  only  derivatives 
of  order  less  than  or  equal  to  two  (case  j  _  -1).  nevertheless,  our  proof  applies  as  well 
to  the  case  j  >  -1,  without  further  difficulties  or  technical  calculations,  as  shown  in 
the  last  part  of  section  3;  the  details  will  be  given  in  a  forthcoming  note,  by 
A.  Def ranceschi .  To  be  more  specific,  here  we  will  concentrate  our  attention  on  the  case 
j  ”  -1 ,  since  our  proof  turns  out  to  be  (slightly)  more  complicated  just  for  this  case. 

Section  4  is  concerned  with  barotropic  motions,  described  by  system  (4.1).  The  main 
goal  of  this  section  is  the  study  of  the  incompressible  limit  of  a  family  of  compressible 
barotropic  flows.  We  assume  that  the  state  function  Px(P)  depends  on  a  parameter  A, 
such  that  the  Mach  number  becomes  small  as  \  goes  to  +”• .  We  refer  to  [5],  [6]  for  the 
physical  justification  of  our  assumptions.  We  will  prove  that  the  solution  of  the 
incompressible  Wavier -Stokes  equation  (4.15)  is  the  limit  of  the  solutions  of  the 
compressible  Navier-Stokes  system  (4.14),  as  A  ♦  +*».  A  similar  result  was  proved  for 
n  <  3,  j  -  1,  p  -  2,  in  reference  [2 j .  The  proof  given  here  applies,  without  supplemental 
difficulties,  to  the  case  W^,p,  if  j+2  >  n/p,  T  £  C3+l,  p(p )  £  C3+i.  The  details  will 
be  given  in  the  forthcoming  note,  referred  to  above. 

finally,  in  section  5  we  will  study  the  existence  of  equilibrium  solutions  (solutions 
such  that  u(x)  •  0,  V  x  £  fl),  for  arbitrarily  large  external  forces  f  »  VF.  In 
particular,  we  will  show  that  if  p(p)  »  RpY,  R  >  0,  T  >  0,  Y  j*  1,  the  equilibrium 
solution  does  not  exist  in  general  (even  for  a  constant  external  force  f).  If  y  *  1, 


-3- 


tnen  the  equilibrium  solution  exists  for  every  bounded  potential  F(x). 

For  small  external  forces,  it  was  proved  in  [11]  that  the  stationary  solutions  are 
stable  (for  the  equilibrium  solutions  this  was  proved  in  [7])*  It  would  be  interesting  to 
study  the  stability  of  the  equilibrium  solutions  for  large  external  forces,  even  if  n  •  1. 
It  is  worth  noting  that  the  core  of  this  paper  is  the  study  of  the  linear  system 

(2.1) ;  see  theorems  2.1  and  3.3.  In  this  system  the  unknowns  are  u  and  o,  the  vector 

field  v  being  fixed.  Due  to  the  term  v*Vo,  system  (2.1)  is  not  an  elliptic  system  in 
the  sense  of  Aqmon,  Douglis,  and  Virenberq  [1;  part  II],  except  if  v(x)  vanishes 
identically  in  fl.  In  fact  (assume,  for  simplicity,  that  u  =*  k  =■  m  «  1,  and  v  <*  0), 

consider  in  SI  the  system  -Au  +  Va  «  F,  div  u  +  v*7 a  +  a(x)a  -  g.  By  using  the  notations 

of  [1;  part  II]  one  has;  L(x,C)  -  vWClCl2"  if  v(x>  *  0;  L(x,C  )  -  ( 1+a(x) )  |5  | 2,1  if 
v(x)  «  0.  Consequently,  the  ellipticity  condition  (see  [1],  eq.  (1.5))  "L(x,C)  t  0  for 
real  £  M  0" ,  cannot  be  satisfied  unless  v(x)  »  0,  for  all  x  e  SI. 

If  v  «■  0  and  if  a  is  small  enough  (for  instance,  if  |a|„  <  1),  the  system  is 
elliptic.  However,  this  last  property  cannot  be  used  to  treat  the  term  vTo  as  a 

perturbation  term.  In  fact,  by  assuming  that  a  c  W1  'p,  the  term  v7o  belongs  just  to 

iP.  In  this  situation,  equation  div  u  +  v*7 a  +  a(x)c  «  g  yields  div  u  e  iP,  and 
equation  -Au  +  Vo  »  F  gives  u  e  W1,p,  o  €  I P.  Hence,  from  the  point  of  view  of 
regularity,  we  lose  one  derivative. 

We  also  point  out  that  the  evolution  problem  is  easier  to  solve  due  to  the  presence  of 
the  term  (3o/3t)  +  v*Vo  in  the  evolution  counterpart  of  equation  (2.1>2. 

Imt  us  briefly  explain  the  main  ideas  utilised  in  the  sequel  (see  also  [2])  to  solve 
the  linear  system  (2.1):  By  applying  the  divergence  operator  to  both  sides  of  equation 

(2.1) 1  one  shows  that  <w+v)A  div  u  -  kAa  -  div  F,  and  by  applying  the  Laplace  operator  to 
both  sides  of  equation  (2.1)j  one  gets 

(1.7)  2L.  4a  +  V'VA 0  -  G(F,g,C)  , 

where  G(F,g,o)  is  given  by  (2.5).  In  order  to  be  able  to  solve  equation  (1.7)  for  Ao  , 

-  1  #P 

we  replace  G(F,g,o)  by  G(F,g,T),  where  t  c  W  is  an  arbitrary  function.  This 

yields  equation  (2.4),  where  (heuristically )  *  should  be  regarded  as  Ac  (actually 


X  «  Ac ,  if  t  «  a  is  a  solution  of  (2.1)).  Once  X  is  known,  we  want  to  determine  div 
u.  By  applying  the  divergence  operator  to  both  si*es  of  (2.1)^,  we  get  equation  (2.10) , , 
where  (heuristically)  X  should  be  regarded  as  div  u  (we  will  show  that  6  -  div  u,  in 
the  event  that  we  have  a  fixed  point  a  -  t ) .  The  remarkable  property  div  u  -  0  on 
r ( 1 )  suggest  us  to  impose  on  9  the  boundary  condition  (2. 10)j.  Once  8  is  known, 
equation  (2.1 )^  gives  u  and  o,  by  solving  the  Stokes  linear  problem  (2.13).  In 
equation  (2.13)  w«  replace  8(x)  by  9q(x)  •  9(x)  -  6,  where  9  is  the  mean  value  of 
8(x),  since  the  compatibility  condition  9^  »  0  is  required  here. 

For  each  fixed  pair  (F,g)  the  above  sequence  of  maps  define  a  map  t  -*■  a  .  We  will 
prove  that  this  map  is  a  contraction  (in  a  suitable  ball),  and  that  the  pair  (u,a), 
corresponding  to  the  fixed  point  a  •  t,  is  a  solution  of  (2.11). 

We  end  this  section  by  calling  the  reader's  attention  to  the  following  papers,  in 
which  results  directly  related  with  those  of  section  3  can  be  found:  M.  Padula  [8] , 

A.  Valll  (11),  A.  Valli  and  W.  Zajaczkowski  (13),  H.  Beirao  da  Veiga  [2],  A.  Valli  (12). 

The  corresponding  evolution  problem  has  been  studied  by  many  authors.  Here,  we 
mention  only  a  sequence  of  papers  by  Mataumura  and  Nishida  (see  [7],  and  references)  and 
the  interesting  Valli 's  paper  (ffj. 

To  readers  interested  in  the  incompressible  limit  of  compressible  fluids,  we  suggest 
papers  by  klainerman  and  Wajda  (5] ,  Wajda  (6) ,  and  Schochet  [9] . 


(*)  Equation  (1.5)j  Implies  that  every  stationary  solution  u  must  verify  the  equation 
div  u  -  0  on  T . 


2.  *Bie  linearised  system 


In  this  section  we  study  the  linear  system 

-piu  -  vV  div  u  +  kVa  «  F(x)  , 

a  div  u  +  v*7'*  +  o  div  v  ■  g,  in  Q  , 


by  adapting  to  the  iP-case  the  method  introduced  in  [2].  We  start  with  some  notation.  We 

denote  by  W^ ,p,  j  an  integer,  1  <  p  <  +“,  the  Sobolev  space  W^'Pffl),  endowed  with  the 

usual  norm  I  ,  and  by  |  |p,  1  <  p  <  +■,  the  usual  norm  in  iP  *  lP(S)).  Hence, 

I  I.  “|  1_.  For  convenience,  we  also  use  the  symbol  W^'P  to  denote  the  space  of 
0#P  “ 

vector  fields  v  in  fl  such  that  v^  e  W^'P(fl),  i  »  1,2,...,n.  This  convention  applies 

to  all  the  functional  spaces  and  norms  utilized  here.  For  j  >  1  we  define 

wJ-P  »  {v  «  W^'P  :  v  «  0  on  T} .  Note  that  wj'P  »  ,p  O  Wg'P,  is  not  the  closure  of 

£><())  -  C^(f»).  Furthermore,  we  set  W  i,p  -  {t  e  Wj'P  :  T  «  0}  ,  VQi,P  -  W^'P  O  W  3,p, 

j  >  1,  where  in  general  J  denotes  the  mean  value  of  A(x)  in  0.  Finally,  for  vector 

fields,  we  define  wJ'P  -  {v  e  wj'p  :  div  v  -  0  on  T) ,  j  >  2. 

0 ,  a  0 

With  the  only  exception  for  section  5,  c,  c^,  i  >  0,  will  denote  positive  constants 
depending  at  most  on  0,  n,  p.  The  symbol  c  may  be  utilized  (even  in  the  same  equation) 
to  Indicate  distinct  constants.  In  section  5,  c  will  denote  an  integration  constant. 

We  denote  by  q  -  p/(p-1)  the  dual  exponent  of  p,  and  by  r*  the  Sobolev  embedding 
exponent  r*  “  (n-r)/nr.  Recall  that  W^'P  ♦  L* ,  if  p  >  n,  and  w1,p  ♦  iP* ,  if 
1  <  p  <  n. 

Now  we  state  the  main  result  in  this  section  (see  also  theorem  3.3): 


THEOREM  2. 1. 


constants  c  and 


assumption 


Let  p  >  n  be  fixed,  and  let  F  e  lP,  g  e  W01,p.  There  exist  positive 
Y,  Y  defined  by  equation  (2.22),  such  that  if  v  €  W^'P  verifies  the 


(2.2) 

then  there  exists  a  unique  solution 

(2.3)  lilul,  +  klol 

2,p  i,p 


Ivl,  „  <  Y'c 
2,P 

« "J:, " 5  ,,p 

‘"(’♦“ife1)  i-v 


of  problem  (2.1). 

B  nihil,*.  . 

m  1,P 


Moreover , 


-6- 


\m 


m 


Proof .  Let  t  e  W  "K,  and  consider  the  linear  problem 


X  +  v*7X  -  G(F,g,T  )  , 


where  by  definition 


(2.5)  G(F,g,T)  -  Ag  +  j  div  F  -  [27v  :  72t  +  Av*7t  +  A  (t  div  v)  ] 

2  2 

Here.  7v  t  7  x  ■  EfSv^/Sx^ )  (3  t/Sx^Sx^).  Note  that  G  is  a  linear  map  from 
Lp  x  wo1,p  x  w  1,p  into  W-1,p  (the  dual  space  of  wj'q),  and  that 

(2.6)  IGI  .  <  Igl.  ♦  -2-  |f|  +  clvl  ItI 

-1,p  1,p  u+v  'p  2,p  1,p 

In  [3]  we  prove  that  there  exists  a  positive  constant  such  that  if 

Ivl,  <  c.  mk/(u+v),  then  there  exists  a  linear  continuous  map  L  :  W-1,p  +  w“1,p  sue 

2  #  p  i 

that  X  «  lg  is  a  weak  solution  of  equation  (2.4),  for  every  G  e  w“1,p.  Moreover 

(2.7)  I  .  <  eld  , 

y+v  -1,p  -1,p 

By  a  weak  solution  of  (2.4)  (see  (3]  for  details),  we  mean  a  distribution  X  e  w-1. 


such  that 


<— ♦  -  div(ev),  X  >  -  «p,G>,  Wee  (Cl) 


Here,  <•«•>  denotes  the  duality  pairing  between  Wq'^  and  W-1'p.  The  above  result 


yields 


Now  let  0  e  W«,p  be  the  solution  of  the  Dlrichlet  problem 


(2.10) 


Clearly, 


(2.11) 


Next,  define 


(2.12) 


(y+v)A0  »  kX  -  div  F  in  Cl  , 


e,r  -o  . 


(utv)l0l  <  c|f|  +  c^(lgl,  +  Ivt,  III,  ) 
1,p  ’p  m  ’  1,p  2,p  1,p 


0 0 < x >  -  S(x)  -  9  . 


Since  0q  •  0,  there  exists  a  unique  solution  (u,a)  e  Wp'P  x  h  1,p 


of  the  linear 


Stokes  problem 


(2.13) 


-uAu  +  k7o  -Ft  v70 


div  u  »  Bg,  in  Cl  , 


U|.  -  0  . 


M 


BS 


Moreover  (see  [1],  (4]) 


(2.14) 


)  u,u,2,p  +  k,0l1,p  *  c(!Ptp  +  <u  +  lv  l.p'  • 

Since  (2.2)  holds,  and  1 0 _ t  ^  <  cl8l.  ,  one  easily  gets 

0  1,p  1,p  1  * 


(2.15) 


lilul  +  klal  +  c(  1  +  u-*l.v-l)  |  f|  +  c  —  ^  I  gt 

2,p  1,p  2  1,p  u+v  p  2  m  ’  1,p 


At  this  point,  we  call  attention  to  the  sequence  of  linear  maps 


<F,g,x)  +  <F,G)  +  (F,X)  +  (f,8)  ♦  (F,80)  ♦  (u,o)  , 

where  F  is  left  unchanged,  and  the  elements  G,  X,  6,  8q,  (u,o),  are  defined  by  equations 
(2.5),  (2.4),  (2.10),  (2.12),  and  (2.13),  respectively.  The  product  map  is  linear  and 


continuous,  by  (2.15).  Hence,  if  (u,,o,)  is  the  solution  corresponding  to  the  data 


( F,g,T , ) ,  it  follows  that  (u-u,,o-o,)  is  the  solution  corresponding  to  the  data 
(0,0,  T-T.):  Consequently ,  10-0,1,  <  ( 1/2  )It -r  ,1  ,  ,  and  the  map  T  ♦  o  is  a 

contraction  in  w  1,p.  Hence,  it  has  a  (unique)  fixed  point  o  »  t. 

Now,  we  prove  that  the  pair  (u,c),  corresponding  to  the  fixed  point  a  «  T  is  a 
solution  of  (2.1).  The  main  point  is  to  prove  (2.1)2,  since  (2.1),  and  (2.1)3  follow 
immediately  from  (2.13).  From  (2.10),  we  get 
(  2. 16)  X  ■  A  div  u  +  div  F  , 

since  A0  »  A0g  =■  A  div  u.  On  the  other  hand,  by  applying  the  divergence  operator  to  both 
sides  of  equation  (2.13),,  and  by  using  (2.13)2  and  (2.16),  we  show  that  X  *  Ac  .  Now,  we 
replace  X  by  the  right  hand  side  of  (2.16)  in  the  first  term  on  the  left  hand  side  of 
(2.8),  and  by  Act  in  the  second  term  on  the  left  hand  side  of  (2.8).  This  yields 
<<p,mA  div  u>  -  <div(gv),ACT>  » 

(2.17) 

=*  <e,Ag  -  2Vv  :  V^a  -  Av»7ct  -  A  (o  div  v)>,  V?  €  V  (0 )  . 

We  claim  that 

2 

-  <div(ev)  ,  Act>  +  <g,27v  :  7  ct  +  Av»7ct>  « 

(  2  e  1  8  ) 

»  <?,A(v*7o)>,  Vy  €  P(0)  , 

where  the  pairing  on  the  right  hand  side  denotes  the  duality  between  P(fl)  and  V  1  (!) )  • 

Let  CTn  c  C3(0) ,  on  ♦  o  in  W1,p(fl),  as  n  ♦  .  The  identity  (2.18),  for  the  functions 

2 

a_,  follows  from  the  formulae  A(v«7o  )  -  v«7Ao  +  27v  :  7  o  +  Av*7o  .  By  passing  to  the 
n  n  n  n  n 

limit  as  n  *  +•,  we  prove  (2.18)  for  o'. 


Proa  (2.17)  and  (2.18)  we  get 


(2.19)  A(m  dlv  u  +  m0  +  v*7 a  +  o  div  v  -  g)  «  0  , 

In  0.  Meuse  that  one  has 

(2.20)  m  div  u  +  ra9  +  v*7o  +  a  div  v  -  g  -  0 

in  fl.  Then,  by  integrating  in  0  both  sides  of  this  equation,  one  easily  concludes  that 
6-0.  Hence  ( 2 . 1 ) 2  holds,  as  desired.  Let  us  prove  (2.20)(  )•  Set 

1)(  ■  {i  <  11  :  dist(c,D  >  c)  ,  -  {*  e  fl  :  dist(s,r)  -  e} ,  where  e  >  0  is  assumed  to  be 

enough  small  so  that  is  a  regular  manifold.  Let  ge(x,y)  be  the  Green's  function  for 

the  Dirichlet  problem  in  fl£,  and  let  x  e  Q  be  fixed.  We  want  to  prove  (2.20)  at  x. 

For  e  e  ]0,(1/2)dist(x,r)(  one  has  sup|Vgc (x,y) |  <  K,  for  y  e  r£  ,  where  K  does  not 
depend  on  e.  Let  0  denote  the  left  hand  side  of  (2.20).  Since  U  is  harmonic  in 
0e,  one  has  U(x)  -  J  Oge(x,y)/3ny)  U(y)dTe(y).  This  shows  that  |u(x)|  <  K/|u(y)|<Te. 
Consequently,  if  we  show  that  lim  inf  /|u(y)|dr£  -  0,  as  E  ♦  0,  it  must  be  that 
U(x)  -  0.  Moreover,  it  is  sufficient  to  verify  that 

(2.21)  lim  inf  ess  /P  |v7c|ar_.  ■  0  , 

e*0  rE  E 

Since  the  function  m  div  u  +  m0  +  a  div  v  -  g  e  wj,p(0)  converges  uniformly  to  0  as 
dlst(y,T)  ♦  0. 

Assume  by  contradiction,  that  there  exist  positive  constants  e Q  and  6  such  that 
the  integral  on  the  left  hand  side  of  (2.21)  is  greater  than  S,  for  a. a.  e  e  10,e0l. 
Denoting  by  K0  the  Lipschitx  constant  of  v  in  0 ,  and  recalling  that  v  vanishes  on 
T,  it  follows  that  /  |Vc|drE  >  «/(KnE),  for  almost  all  e  e  ]0,e0(.  This  contradicts 
the  estimate 

e 

/0  </r  |Va|dTE)dE  -  Jn/n  |Va  |dy  <  d-  . 

E  E  A 


(2)  In  case  that  j  >  0  (see  theorem  1.1)  the  proof  of  (2.20)  is  immediate,  since  the  left 
hand  side  of  (2.20)  is  harmonic  in  0  and  vanishes  on  T  (since  dlv  u  -  ,  v-  0, 

div  v-g-0,  on  T).  However,  if  j  -  -1,  this  argument  has  to  be  carefully  handled,  since 
the  function  7c  has  not  a  trace  on  T . 
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The  existence  pert  of  theorem  2.1  is  completely  proved.  As  in  (2),  section  2,  we  show 
that  there  exists  a  constant  Cg  such  that  the  solution  of  (2.1)  is  unique  if 
Ivl  <  (mu  k)/(c  (m  +  !v  I  )  21  ,  where  u.  =  min{u  ,g  +v}  . 

Z/p  U  0  0 

For  convenience,  we  set 


(2.22) 


Cim  m  o 

2c2(™+,V,)'  c0(U  +  M)2' 


Then,  inequality  (2.2)  includes  all  the  assumptions  on 

2.  1. 


Ivl  utilized  in  proving  theorem 

2<P 
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In  this  section  we  prove  the  existence  and  uniqueness  theorem  1.1  in  case  that 


j  «  -1,  l.e.  we  prove  theorem  3.1  below.  We  close  the  Section  by  showing  how  to  adapt  the 
proofs,  in  case  the  j  >  0. 

Recalling  (1.4),  we  set 


<^(0,0)  -  (tf  (T  ,B)  I 


T  -  sup  !«I  (<?,<*>  I  ,  S  -  sup  — nr—; - y—. rr 

1  Id,!,)  1  1  1(1,1,)  ,(°'a)  (T'B> 


,1-1,2  . 


Moreover,  we  set  u>,(t,s)  -  sup  |ii) ,  (<J  ,a  )  |  ,  for  |o|  <  t,  |a  |  <  s.  We  denote  by  c3  a 
positive  constant  such  that 


(3.1) 


T  _  <  C,lTl 


1#p  '  I®!-  *  c3,8,1,p  '  **  e  M  1,P'  «  «  w 


For  convenience.  In  this  section  we  denote  by  c’,  c^*,  1  >  0,  positive  constants 
depending  at  most  on  (1,  n,  p,  v,  \>,  h,  m,  cQ,  c v,  x»  XQ.  X,.  1*  t,»  Tit  Sj_.  The  symbol 
c*  may  be  utilised,  even  in  the  same  equation,  to  denote  distinct  constants. 

One  has  the  following  result)3): 


Theorem  3.1.  Let  p  >  n  be  fixed,  there  exist  constants  c0,  c, ,  such  that  if 
(  t  if,  g  <  w  1,p,  h  t  iP,  and 

(3.2)  If |p  +  '9*t#p  ♦  lhlp  <  c'0  , 

there  exists  a  unique  solution  (u,o,a)  e  W2,p  x  W  1,p  x  Wg,p  of  problem  (1.5),  in  the 
hall 

(3.3)  lu«_  e'  • 

lip  1 

Hence,  there  exists  a  unique  solution  (u,p,c)  c  w2,p  *  w1,p  *  w2,p  of  problem  (1.1), 

(1.3) ,  in  the  corresponding  neighborhood  of  (0,m,^Q). 

In  order  to  prove  theorem  3.1  we  state  some  auxiliary  results.  The  verification  of 
the  following  statement  is  immediate,  and  is  left  to  the  reader: 

Lemma  3.1.  Let  t  «  W  1,p,  B  e  wj,p,  verify  the  assumptions 


(3)We  assume  here  that  p(P,C)  has  Lipschitr  continuous  first  order  derivatives.  However, 
theorem  3.1  still  holds  if  p(p,t)  is  only  assumed  to  be  continuously  differentiable. 


',1,1 


m 


and  let  v  e  W2,p,  f  c  iP ,  g  c  W0',P»  h  «  I?.  Then 


u 

SI 


|F(f,v,T,8)|p  <  f  »(!f|p  +  o.vl^p)  +S1(|tL,|8L).tI1#p  ♦  T2»8I1>f 


*V.il 


(3.5)  <  H(h,v,t,8)  <  =r  a  h  +cc  ml  vi  .  181.  + 

1  'p  2  '  'p  v  1,p  1,p 


l+  c(x0+x1),v,1,p,v,2,p  +  c  —  (,V,2/P,T,UP  +  ,g,P)  * 

Theorem  2.1,  lemma  3.1,  and  claaslcal  results  for  the  DLrlchlet  problem  (3.7),  yield 
the  following  statement: 

2  p 

Theorem  3.2.  Let  v  €  Wjj'p  verify  (2.2),  and  assume  that  the  hypotheses  in  lemma  3.1 
hold.  Then,  there  exists  a  unique  solution  (u,o)  e  wg,P  *  W  of  problem 

-liAu  -  vV  div  u  +  kVa  -  F(f,v,T,B)  , 

<  m  div  u  +  v*7c  +  a  div  v  «  g,  in  n  , 

(3*6)  — 

utr  -  0  , 


and  a  unique  solution  a  e  Wg'P  of  the  problem 


Moreover , 


-xAa  »  H(h  ,v,x  ,8  ) ,  in  0  , 


a,r  -  0 


u.ul2  p  ♦  k.a.,#p  <  c(1  ♦  HJM,  i«|f|p  ♦  lv.J  p  ♦ 

♦51(M..I«lt>).ti1#p  +  T2.e.1fpi  , 

X'a,2,p  ‘  C,h,p  *  °  Cv",V,1,p'8,1,p  +  C<X04X1,,V,l.p  + 


+  c  — —  (Ivl,  It  I ,  +  g  ) 

m  2,p  1,p  '^'p 


88L 


in  the  sequel,  the  solution  (u,o,a)  of  problems  (3.6),  (3.7)  with  data  (v,t  ,8  )  is 
denoted  by  (u,o,a)  -  T(v,t,8).  Let  us  write  (3.8),  (3.9)  In  the  abbreviated  form 


I,  _  ♦  „  <  ci(|f|  ♦  Igl  ,  +  Ivl2 

2»P  >,p  2  P  1,p  2,p 


(3.10) 


'<*•,  n  ‘  C2l'h,D  +  r.  +  „  +  »tl*  „  +  181*  ) 

2#  P  2  p  i#p  2,p  1#p  1,p 


(3.11) 


-  mln{yk,  ,  (?  )  >  25  c,  ,  5  cid+cl^ 


let  r  e  ]0,rn),  and  assume  that 


(3.12) 


and  that 


(3.13) 


l*lp  <  *2,  Il»!p  <  r2,  Igt,  p  <  r2  , 
lvl_  <  r,  Itl  <  r,  I8»,  <  5  cl  r2  . 

2,p  1,p  1,p  2 


Wider  these  assumptions,  the  conditions  (2.2),  (3.4)  are  fulfilled.  Consequently, 

from  (3.10)  it  follows  that  (u,0,a)  "  T(v,t,B)  verifies 

(3.14)  lul,  <  r ,  lal  <  r,  lal  <5  c*  r2  . 

2,P  1,p  1,p  2 

Consequently  T(Br)  C  Br,  where  Br  is  the  subset  of  W^,p  *  W  1,p  *  w2'P  defined 

by  equations  (3.13).  In  order  to  accomplish  the  proof  of  theorem  3.1,  we  will  show  that 

T  is  a  contraction  in  Br ,  respect  to  the  norm  Wg*2  *  L2  *  Wg'2,  if  r  is  sufficiently 

small,  it  follows  then,  by  the  contracting  mapping  principle,  that  T  has  a  (unique) 

fixed  point  (u,o,a)  “  (v,t,8)  in  Bf  (we  could  also  use  the  Schauder's  fixed  point 

12  2  12 

theorem,  since  8r  is  a  convex  and  compact  subset  in  respect  to  the  x  L  x  w  ' 

topology) . 

Let  (a, a, a)  -  T(v,r,8),  (u1,o1,o1)  «  T(  v^  ,t  ^,8  , ) ,  F  -  F(f,v,t,8),  F,  ■ 


F( f  ,v ^  ,t  ^ ,8 ^ ) ,  H  *  H(h,v,T  ,8) ,  ■  H(h,v^,x  ^  ,8  1 ) .  One  has,  in  0, 
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(3.15) 


(3.16) 


-uAfu-u^)  -  vV  divtu-u^)  ■*•  kVfa-o^)  -  T  -  , 


m  div(u-u^)  +  v^*V(o-«^)  +  (v-v^Va  + 


+  divfv-v^)  +  (cr^)  div  v  «  0  , 


-X  Afa-a^  «  H  -  H,  . 


Here,  we  use  the  notation  1  1^  j  «  I  1^.  By  multiplying  both  aides  of  equation 
(3.15),  by  m(u-u,),  both  sides  of  equation  (3.15)2  by  k(o-o,),  by  integrating  the 
resulting  equation  in  ft,  and  by  adding  them  side  by  side,  we  show  that 


(3. 17) 


m(y-|  v  | )  IVtu-u^lg  <  mlP-P^I  lu-u^l  + 


♦  cklo.l.  l»-v,l  .lo-o.l.  +  cklvl.  I9-O.I. 

1  1,p  11  10  2,p  1  0 


In  proving  (3.17)  we  utilize  some  well  known  Sobolev  inequalities.  Note  that  if 
2*  -  2n/(n-2) ,  one  has  (1/2*)  +  (1/2)  +  (1/p)  <  1. 

Arguing  as  on  proving  equation  (3.12)  in  reference  [2),  we  show  that 


(3.18) 


Moreover , 


(3.19) 


lu-u  I2  +  c*(1  -  c’lvl  )  la-o  |2  < 
11  4  5  2,p  10 


<  c' la  I  2  I v-v  I  2  +  c' I  M  l  2  . 

1  1,p  11  1-1 


la-a,!  1  <  c'IH-H1l_1  . 


Straightforward  calculations  show  that  (see  appendix) 


r-V.,  <c*(lf|  *ivi.*  )It-«i,o* 


(3.20) 


+  c’(1  +  ,Tl1,p><'v,1,p  +  '»l'l,p"”,',  * 

*  C’(,T,1.P  ♦  ,Ti'i,p  *  '8,1,p  ♦  ,8l'l,p>(,T-Tl'o  *  ,Ml'o>  ' 


and  that 


f 

tw_ 

iuV 


(3.21) 


IH-H.I  ,  <  C*(lvl,  _  +  Iv.l,  ♦  III,  +  161  .  )  lv-v,l,  + 

1  “ I  '  #p  '  •  #p  '  *p  l»p  1  I 

+  c-(|h|p  +  Ulp  ♦  C8«1>p  +  mi>p,.v.1jp  ♦  (1  ♦  .T,.1#p  *  I6,.1>p  + 
+  ,Ti'?.P  *  ,6i,i,P),vi,i.p1,T-ti,o  +  c,('^p  +  ,v,1.p'T,1,p  + 


►  u 

I 


+  'Vi,p,Ti'i.p  +  'Vi^'8"6^  •  • 

By  using  (3.12),  (3.13)  it  follows  that  the  coefficients  of  Iv-v^  Ix-^l., 

IB-6  ^1 g,  in  the  right  hand  sides  of  equations  (3.20),  (3.21),  are  polynomials  on  r, 

vanishing  for  r  «  0,  and  with  coefficients  which  are  positive  constants  of  type  c' . 

2 

Similarly,  the  coefficients  c'lvl  and  c'lo  I  (that  appear  in  equation  (3. IS))  are 

5  2,p  i  i,p 

polynomials  on  r  of  the  above  type.  Since  the  exact  form  of  these  polynomials  is  not 
important  here,  we  will  denote  them  by  the  symbol  e  -  e(r). 

By  using  the  above  notation,  the  estimates  (3.18),  (3.19),  (3.20),  (3.21),  yield 


(3.22) 

(3.23) 
Hence, 

(3.24) 


lu-u1l^  (c'-ejlcr-cr^l  *  <  elv-v^l*  ♦  elx-x^l*  +  elB-B^g 
la-a^l^  <  elv-v^l*  ♦  elx-x^l*  +  elB-B^i* 


lu-u^*  ♦  (c'-eMo-c^*  ♦  Ite^l  *  < 


<  elv-v^  ♦  e.x-x,^  *  e,B-ei'? 


Hence,  for  a  sufficiently  small  value  of  r  (depending  only  on  0 ,n,p, . . . ,T^,SL) ,  T  is  a 

12  2  12 

contraction  in  Br,  with  respect  to  a  suitable  norm  in  wq'  *  L  *  Mo'  •  Th*  proof  of 
theorem  3. 1  is  accomplished. 

□ 

Proof  of  theorem  1.1.  The  changes  to  be  made  on  the  proof  of  theorem  3.1,  in  order  to 
adapt  it  for  all  values  of  j,  are  quite  obvious.  He  start  by  showing  it  for  theorem  2.1, 
which  is  the  main  tool  in  our  proofs. 

The  proof  of  theorem  2.1  was  done  by  solving  the  sequence  of  linear  problems  (2.5), 
(2.4),  (2.10),  (2.12),  (2.13),  and  by  getting  suitable  estimates  for  the  corresponding 
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%  > 


solutions .  Let  us  show  that  the  proof  goes  on,  in  case  that 
j  -  -t. 


j  >  0,  exactly  as  for 


Problem  (2.5):  The  following  estimate  for  the  solution  of  problem  (2.5)  is  easily 

obtained  by  using  Sobolev's  inequalities: 

IGI  <  Igl.,  „  +—  IFtK.  +  Clvlj  ,  ITI^, 

j,p  3+2, p  U+V  j+1,p  3+3, P  3+2. p 

Problem  (2.4):  The  following  result  is  proved  in  [3]:  There  exist  positive  constants 


c  and  c.  (here  the  constants  depend  also  on  j)  such  that  if  Ivl 


<  c  mk/(u  +v ), 


- - - "  — - - * -  J'  . . 3+3»p  '  “1 

and  v  «  0  on  V,  then  there  exists  a  solution  X  of  problem  (2.4),  moreover  one  has 

— 111  .  <  eld  . 

2(y +v)  j,p  j,p 

Problem  (2.10):  Por  the  solution  8  of  the  Dirichlet  problem  (2.10)  the  estimate 

(U+v )  191  <  c(IPI  +  kill  ) 

3+2, P  3+>#P  3»P 

is  well  known. 

Problem  (2.12):  Obviously,  the  function  8q,  defined  by  equation  (2.12),  verifies 

the  estimate  < 9 _ I  .  _  <  clfll ... 

0  j+2,p  3+2.P 

Problem  (2.13):  The  solution  (u,o)  of  this  Stokes  linear  problem  (see  [1],  [4]) 
verifies  the  estimate 

U,U,j+3.p  +  k,<,'j+2,p  *  C<,Mj+1,p  +  (lJ*lVl),90'l+2.p>  * 

By  putting  together  the  above  estimates,  one  easily  gets  (c.f.r.  (2.15)) 

U  lul  ,  +■  klol  J  ,  <  ^IH.  ,  ♦ 

j+3,p  3+2. P  2  3+2. P 

+  C(1  +HlbLL)  +c  >L!MIq,  . 

U+v  j  +  1,P  2  m  3+2, P 

By  arguinq  as  in  section  2  (recall  the  footnote  (1))  one  proves  the  following  result: 

Proposition  3.3.  Let  p,  j,  and  T  be  as  in  theorem  1.1  and  let  Pc  W^+1'p, 
g  e  ^'P.  Then,  there  exist  positive  constants  c  ”  c(fi,n,p,j),  Y  “Y(fl,n,p,3»U»v,m), 


such 


that  if  ve  wj+2'p  verifies  the  condition  Ivl ,  .  <  Yk,  then  there  exists  a 

-  o,a  —  ■■  ■  . . . . .  3+3, p  ■■  . — — 


unique  solution  (u,o)  c  W^+*'p  *  W^+2'P  of  problem  (2.1).  Moreover , 

ulul,  .  +  klol  .  „  <  c<1  +  U.*.LVI)IPI 

j+3,p  3+2, p  U+v  3+1, P 


+  c  V*\?±  ,  g, 


3+2, p  • 
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The  rest  of  th«  proof  of  Theorem  1.1  (th«  "non-lin ear  part*)  cloaaly  follows  the  proof 
of  Theorem  3.1.  For  an  arbitrary  j  >  -  1,  we  define  the  set  Br  by  the  Inequalities 
(c.f.r.  (3.13)) 

Ivl.  ,  <  r,  ItI .  ,  <  r,  161.^,  <  c'r2  , 

j+3»p  .  j+2,p  j+2,p 

and  we  prove  that  (for  sufficiently  small  values  of  r)  T(Br)  C  Br>  and  that  T  Is  a 

contraction  in  Br  with  respect  to  the  norm  W^'2  *  L2  »  W^'2.  The  proofs  of  these  two 

statements  in  case  that  j  >  -1,  differ  from  the  proofs  given  above  only  by  the  particular 
Sobolev's  inequalities  to  be  used. 

Before  ending  this  section  we  remark  that  in  Theorem  1 . 1 ,  if  h  e  'P  ( instead  of 

h  e  »i+1'p)  then  C  <  wl+2'p,  and  if  h  «  then  C  e  M1+4,p. 
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4.  The  incompressible  limit. 


In  this  section  we  study  the  system 


-liAu  -  v7  div  u  +  7p(p)  -  p  [f  -  (u»7)u) 


div(pu)  “0,  in  ft  , 


u)r  -  o  . 


describing  the  barotropic  motion  of  a  compressible,  viscous  fluid.  Since  we  are  interested 


on  studying  the  limit  of  the  solution  u  when  k  =  p^  (m)  ♦  ,  it  is  necessary  to  state 


an  existence  result  for  prohlem  (4.1)  in  which i  (i)  the  dependence  of  some  suitable 


structural  constants  of  the  state  function  p(p )  in  terms  of  k  is  given;  (ii)  the 


dependence  on  k  of  the  constants  appearing  on  the  estimates,  is  shown.  This  is  the  aim 


of  theorem  4.2  below. 


bet  p(p)  be  continuously  differentiable  in  [m  -  (lA)i  m  +  UAllr  where 


0  <  l/k  <  m/2,  and  write  Pp  in  the  form  Pp(m+o)  “  k  ♦  u(o).  We  assume  that  there 


exists  a  e  ]0,1]  such  that 


|u«j}|  <  g|o|“  ,  vo  e  HA.tAl  , 


and  we  denote  by  *  a  positive  constant  for  Which 


S  <  4k1 


Arguing  as  in  proving  theorem  3.2,  we  show 


Theorem  4. 1 .  Let  the  above  assumptions  on  the  state  function  p(p )  hold ,  let 


v  e  WQ  d  verify  (2.2),  and  let 


Iktl  .  <  — 


1»P  c. 


If  f  e  I P,  there  exists  a  unique  solution  (u  ,o )  e  WQ'*j  x  w  ,p  of  the  problem 
(  -pdu  -  v7  div  u  +  kVc  -  F(f,V,T) 


m  div  u  +  v’Vc  +  o  div  v  «  0,  in  (l  , 


uir  “  0  . 


where  F(f,v,t)  -  (x+m)(f  -  (v*7)v]  -u(t)7t.  Moreover , 


'"'a.p  +  'to'l,p  <  ci,f'p  +  CJ'V,1,p  +  C3  * 


Here  the  positive  constants  cj,  cj,  cj,  depend  at  most  on  ft,  n,  p,  v,  v,  m  (and 
are  distinct  from  the  constants  cj,  cj,  cj,  in  section  3).  Following  section  3,  we 
denote  by  (u,a  )  »  T(v,t)  the  solution  of  system  (4.5).  Define 
(4-7)  “  mingle,  r-  ,  ^  1/a}  , 


rQ  -  min{y)t,  ,  ( -- -1^.) 


and  assusMi  that  f  t  tP  is  fixed,  and  verifies 

(4.8)  |f |  <  r/(3  cj)  . 

One  easily  shows  that  T(Br)  C  Br,  where  Bf  is  defined  by  the  inequalities 

(4.9)  I  vl  _  <  r,  I  krr  I  .  <  r  . 

2,p  1,p 

Moreover,  there  exist  positive  constants  c',  cj,  cj,  depending  at  most  on  0,  n,  p,  u, 
m,  such  that 


(4. 10) 


lu-u^  +  c^(1 
<  c'l  to 


-  c5,v,2,p)'k((,-<,1),0  < 

2  2  2 
,  +  C*IF-F  I  , 

1,p  11  1-1 


where  (u^o,)  - 
inequality  (3.18) 


(4.11) 


T(v1,r1),  F  -  F(f  ,v,T  ),  F1  »  F(f,v1,T1).  Inequality  (4.10)  is  proved  as 
.  Straightforward  calculations  (see  (3.20))  show  that 

i  <  clflp,T-ri'o  +  C<1  +  ,T,i,p)<,v,i.p  +  ,vi,i,p,,v-vi’i  + 

2 

♦  cl  V.l  .  lT-t.l„  + 

1  1,p  1  0 


+  c“^(lkTl“  +  IkT  l“  )lk(T-T  )l  . 

3T  1,p  1  1,p  1  0 

Arguing  as  in  section  3,  one  easily  verifies  that  there  exists  a  positive  constant 

r.,  <  rQ,  depending  only  on  SI,  n,  p,  v,  v,  s,  I,  1,  kg,  such  that  for  r  <  r1#  T  is  a 

12  2 

contraction  in  Br,  with  respect  to  the  Wg ’  *  L  norm  (for  convenience,  we  denote  by 

kg  a  positive  constant  such  that  k  >  kg).  The  above  result  proves  the  following 
theorem: 


Theorem  4.2.  Let  p  >  n  be  fixed,  and  let  the  above  assumptions  on  the  state 
function  p(p  )  hold.  Then,  there  exist  positive  constants  cj,  cj. 


depending  only  on 


n,  n,  p,  p,  v,  m,  l,  ♦  and  kg,  auch  that  If  f  £  1^  verifies  |f|p  <  c£,  then  there 
exists  a  unique  solution  (u,p)  £  Wg'p  *  W1,p  of  problem  (4.1),  (1.3)p  In  the  ball 


lul  «  c*  ,  Ip-ml  <  c*A- 
z,p  7  1,P  7 


From  now  on  we  assume  that  px(p)  Is  a  family  of  state  functions  depending  on  a 


parameter  X  £  (XQ, +<»>[,  and  such  that  kx  ♦  +®  as  X  ♦  +*•.  Our  assumptions  on  (p ) 


are  the  following: 


Let  Dpx(p)  denote  the  derivative  dp^(p)/dp.  We  set  kx  -  Dpx (m) ,  we  assume  that 


kx  >  kg  >  0,  and  we  suppose  that  P\(P)  is  defined  and  continuously  differentiable  in  the 


interval  [m  -  tAx«  *  +  IA\]«  where  0  <  t  <  (kgm)/2.  Moreover,  we  assume  that 


|Dpx(p)  -  Dpx(m)  |  <  Sx|p-m|  ,  for  a  fixed  a  c  ]0,1].  Hence,  by  setting  Dp^  (m-w ) 


kx  +  ux(0),  one  has  ux( 0)  -  0,  and 


(4.12) 


ux  (CT)  |  <  |a  va  £  t-t/k^  ,  t/kx) 


We  assume  that  there  exists  a  positive  constant  4  such  that 


(4.13) 


Sx  <  *  k^**  ,  VX  >  XQ  . 


These  quite  general  assumptions  contain  the  physically  interesting  cases  described  by 


klainerraan  and  Majda  in  (5].  In  particular,  if  k,  a  and  S  are  the  parameters  relative 


to  a  given  state  function  p(p),  and  if  we  define  P\(P)  *  X^Cp),  then  Sj^  <  (S/k)kx  _ 
Hence,  Sx  <  ♦kx+a,  where  ♦  «  S/fX^k '■**).  In  that  example  we  may  view 


M  -  l/X  as  the  Mach  number  (see  [5]).  One  has  the  following  result: 


Theorem  4.3.  Let  p  >  n  be  fixed,  and  let  the  above  assumptions  on  the  family  of 


state  functions  px(P)  hold.  Then,  there  exist  positive  constants  c£,  c$,  depending  at 


most  on  0,  n,  p,  u,  v,  ra,  L,  ♦  ,  kn,  such  that  if  f  £  L^,  | f |  <  ci,  then  the  following 


statements  hold: 


(1)  for  each  X  >  X0,  the  problem 


(4.14) 


-uAuj  -  vV  div  ux  +  Vp^(Px)  -  Px  [f  -  (Uj'V)^)  , 


div(Pxux)  -  0,  in  0  , 


(u, ) I-  -  0, 


aw  i  mr 


ww ww  w i  *jtwj i  mui  i  m\ m' 


’«■'  Wtllliimii 


has  a  unique  solution  (u^,p^)  e  W^'p  *  W1,p  in  the  ball  luj^lj  p  *  cg' 

,p\  -  ■’l.p  4  «yv 

(ii)  _If  lim  -  +*•  as  J  +  than  ♦  q.  weakly  in  Wg'p.  div  ^  *  0  weakly 

In  wj'p,  px  ♦  *  strongly  in  W1,p,  ^Px<P\>  +  7*  weakly  in  tf,  where  (q,»7*)  le  the 
unique  solution  of  tha  incoapreeeible  Havier-Stokes  equations 

-uAu^  +  Vw(x)  -  m[f  -  (u—»V )u—]  , 

div  u  •  0  ,  in  0  , 

(u  )i  -  0  . 

Ir 


(4.15) 


Proof.  Part  (i)  follows  from  theorem  4.2.  Moreover,  the  estiaates  show  that  there 
exists  u.  e  Wg'p  such  that  u^  ♦  u«,  and  px  ♦  «,  in  the  topologies  indicated  in  the 
statesMnt.  the  convergence  of  all  the  sequence  \iy  to  q.  follows  from  the  uniqueness  of 
the  solution  of  (4.15),  which  holds  if  c£  is  sufficiently  small. 

From  equation  (4.14)2,  one  easily  verifies  that  div  ux  ♦  0  in  tP .  Since 
Idiv  u,1  ,  <  c*.  one  gets  the  weak  convergence  in  wj'p.  Clearly,  div  q,  -  0. 

A  l,p  9 

Finally,  we  pass  to  the  limit  in  equation  (4.14),.  One  hae  Vpx(Px)  *  mAH.  + 
m(f  -  (u-*V)u-]  ,  weakly  in  iP,  since:  Au^  ♦  Aq,  and  -vV  div  u^  ♦  0  weakly  in 
tP»  Px  ♦  m  strongly  in  W1,p;  and  PX<Ux*?)ux  ♦  sUq/Vjq.,  strongly  in  I P.  Obviously 
the  limit  of  the  sequence  Vpx(px) 


for  »  €  W  1,P. 


must  be  of  the  form  V*(x) 


5.  Tt>e  equilibrium  solutions 


The  results  described  in  this  section  are  obtained  by  using  very  elementary  methods. 


The  proofs  are  independent  of  the  results  of  the  proceeding  sections. 


By  an  equilibrium  solutions  (e.s.)  of  system  (4.1)  [resp.  (1.1)1  we  mean  a  solution 


(u,p)  (resp.  (u.p.Ol  «uch  that  u  -  0  in  0,  and 


0  <  ess  inf  p(x)  ,  ess  sup  p(x)  <  +» 
xtfi  xcft 


If  h  »  g  «  0  the  two  problems  are  equivalent  (set  p(p)  «  p(P.Cg)  in  system  (4.1)). 


In  the  sequel  we  will  consider  the  system  (4.1).  We  call  equilibrium  solutions  the 


functions  p(x),  in  the  class  (5.1),  such  that  (0,p(x))  solves  (4.1).  It  is  easily 


shown  that  a  necessary  condition  for  the  existence  of  the  equilibrium  solution  is  that 


f  *  VF,  for  some  potential  F.  Moreover,  the  equilibrium  solutions  are  the  solutions  of 


the  equation 


Vp(p(x))  -  p(x)Vf(x),  x  r  n 


Here,  we  are  interested  in  the  study  of  the  e.s.  under  the  effect  of  arbitrarily  large 


external  forces  f  -  Vf.  We  will  present  useful  necessary  and  sufficient  conditions  for 


the  existence  of  the  equilibrium  solution  for  an  arbitrary  Ft  L~  (0  ) . 


In  the  sequel,  p  is  a  continuously  differentiable  real  function  defined  on  R+ 
{s  c  R  s  a  >  0),  such  that  p'(s)  >  0,  W  s  c  R+.  We  look  for  equilibrium  solutions 


verifying  (5.1)  and  (1.3)v  for  a  fixed  m  >  0.  We  define 
(5.3)  Ms)  -  /“  t-1p*(t)dt,  V  s  c  R+  . 

Ilk 

We  denote  by  ]a,b(  the  range  of  *,  ]a,b[  •  »(R+).  One  has  <  a  <  0  <  b  <  •♦», 
since  v (m)  -  0.  We  define  *  »  s-1.  Clearly,  4()a,b()  -  R* •  We  set  4(a)  -  0, 


l>t  o  and  F  be  defined  and  measurable  in  fl,  and  let  p  verify  the  assumptions 


(5.1)  and  (1.3),.  Then  p  is  said  to  be  a  weak-equlllbrlum  solution  if  there  exists  a 


real  constant  c  such  that 


ir(p(x))  “  F(x)  +  c,  a.e.  in  0 


Moreover,  p  is  said  to  be  a  strong-equilibrium  solution  if  p  t  W1'1  and  if  (5.2) 


m 


m 


m 


holds  a.e.  In  (l.  A  veak-equillbrium  solution  is  a  strong  e.s.  if  and  only  if  p  t  w1'1. 
Furthermore,  P  £  W1'1  if  and  only  if  Ft  W1'1.  Hence  it  is  sufficient  to  take  in 


account  the  weak  formulation  (5.4).  Note  that  equation  (5.4)  shows  that  p  and  F  have 


the  same  regularity,  if  p  is  sufficiently  smooth.  For  the  sake  of  convenience  we  will 


concentrate  here  on  the  L  case. 


Definition  5. 1 .  Let  F  e  L  .  A  function  p  is  called  an  equilibrium  solution  of 


system  (4.1)  if  p  £  L*  and  if  ( 5.4) ,  (5.1)  and  (1.3),  hold. 


We  set  mg  •  ess  inf  F  in  ft.  Mg  ■  ess  sup  F  in  ft.  One  has  the  following  result: 


Theorem  5.2.  Let  F  e  L  ,  be  given.  There  exists  an  equilibrium  solution  p(x)  if 


and  only  if  there  exists  a  constant 


c  £  ]  a  —  mg,  b  —  Mg  [ 


such  that 


-pq-  /n  *(c  +  F(x))dx  -  ra 


If  such  a  constant  exists  then  the  (unique)  equilibrium  solution  is  given  by 


0  (x)  -  *(C  F(  x)  )  ,  V  X  £  n  . 


Proof .  The  condition  (5.6)  corresponds  to  condition  (1.3),.  Uniqueness  is  shown  by 


noting  that  the  left  hand  side  of  (5.6)  is  a  continuous  and  increasing  function  of  c  in 


the  interval  (5.5).  The  details  are  left  to  the  reader. 


Theorem  5.3.  Under  the  assumptions  of  theorem  5.2,  there  exists  an  equilibrium 


solution  P(x)  if  and  only  if 


mg  <  b  -  Mg  , 


w  *(a  "  "o  +  F<*))dx  <  ra  <  Irff/n  *(b  -  M0  +  P(x))dx  ■ 


In  this  case  the  e.s.  p(x)  is  given  by  (5.7),  where  c  is  the  (unique)  solution  of 


(5.5),  (5.6). 


Proof »  Condition  (5.7)  follows  from  (5*5).  Moreover,  (5*6)  holds  for  some 


c  c  ]a-mgf  b-Mg [  if  and  only  if 

(5.10)  lira  +  -r^-r  /  *(c  +  F(x)  )dx  <  m  , 

c*(a-mQ) 

(5.11)  lim  "f(TT  +  F(x))dx  >  m 

e*(b-Mor  1 

These  two  conditions  are  equivalent  to  (5.1),  by  a  well  known  generalization  of  Beppo 
Levi's  theorem  to  functions  taking  values  in  [0,+*].  Note  in  particular  that  the  integral 

on  the  right  hand  side  of  (5.9)  is  well  defined,  since  0  <  ♦(b-Mg  +  F(x))  <  +»  ,  a.e.  in 

(1. 

a 

Corollary  5.4.  If  a  =  [resp.  b  *=  +*»]  then  the  conditions  (5.8)  and  ( 5.9)  ^  [ reap . 
(5.9)2]  hold,  for  every  F  c  L°* .  In  particular,  if  ]a,  b[*>]  -  «• ,  +  «•  [  the  equilibrium 
solution  exists  for  every  F  e  l“. 

Corollary  5.5.  Let  F  e  L~ .  Then,  the  equilibrium  solution  exists  if 

(5.12)  M0  -  raQ  <  min{-a,b} 

In  particular  this  holds  if  ] F | „  <  (1/2)  min{-a,b). 

Proof ■  If  (5.12)  holds  then  a  -  mQ  +  F(x)  <  0,  and  b  -  Mg  +  F(x)  >  0,  a.e.  in  !). 

Hence  »(a-mg  +  F(x))  <  m  <  *(b-Mg  *  F(x)),  and  assumption  (5.9)  holds. 

□ 

Remark  5.6.  The  lack  of  condition  (5.9)^  [resp.  (5.9)2]  corresponds  to  the  formation 
of  vacuum  [resp.  points  of  infinite  density].  If  a  »  [resp.  b  -  -t~ ]  the  first  [reap, 
second]  phenomena  does  not  occur,  under  the  effect  of  bounded  potentials  F. 

Let  us  now  consider  the  one  dimensional  case  (l  «  ]0,l[,  t  >  0,.  in  the  presence  of 
external  forces  f  t  L^.  The  potential 

F(x)  -  /*  f  ( t )  dt 

is  then  an  absolutely  continuous  function  in  [0,1].  Consequently,  the  equilibrium 
solution  p  (if  exists)  belongs  to  W1'1  and  verifies  (5.2).  The  sufficient  condition 


(5.12)  can  be  written  In  the  equivalent  form 


(5.13) 


1/5  f ( t ) dt |  <  min{-a,b}  ,  V  x,y  e  [0,t]  . 


particular,  |f| j  <  min{-a,b}  ie  alao  a  sufficient  condition  for  the  existence  of 

the  e.s. 

Examples .  As  shown  before,  the  existence  of  the  equilibrium  solution  under  a  given 
potential  F,  strongly  depends  on  the  values  of  a  and  b.  It  is  of  interest  to  consider 
the  classical  case  p(p)  -  SpY,  R  >  0,  y  >  0.  It  follows  immediately  from  our  results  that 
if  Y  >  1  (reap.  Y  <  1],  vacuum  [rasp,  infinite  deneitiee]  may  occur.  On  the  contrary, 
if  Y  ■  1  the  e.s.  exists,  for  every  bounded  potential  F. 

Let  us  give  an  example ,  corresponding  to  the  main  case  y  >  1.  For  convenience,  we 
consider  the  one  dimensional  case  0  -  ]0,1[.  Let  m  »  1,  Y  »  2,  R  -  1/2.  Then  p(p  )  - 
<1/2)p^,  w(p)  "  p-1,  *(8)  ”  Its,  a  ■  - 1 ,  b  ”  +“.  In  the  light  of  corollary  5.4,  we  have 
only  to  check  the  first  condition  (5.9),  which  corresponds  to  the  formation  of  vacuum.  By 
assuming  that  f(x)  ”  8  is  constant,  this  condition  is  |6|  <  2.  Formulae  (5.6)  and  (5.7) 
yield  c  -  -6/2  and  p(x)  -  1  +  8(x  -  1/2),  respectively.  If  B  >  0,  the  forces  act  on 
the  positive  direction,  and  the  point  x  “  0  is  subject  to  the  largest  decompression. 

For  6  •  2  the  vacuum  is  attained  at  the  point  x  ■  0,  since  p(0)  -  0. 
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jggendlx .  For  the  reader's  convenience,  vie  will  prove  here  the  estimate  (3,21).  T 
proof  of  (3.20)  is  easier,  and  is  left  to  the  reader.  Recall  that  under  the  assumptions 


made  in  section  3  one  has  1  t  <  t,  |t,|„  <  t,  |B  |„  <  t,,  <  tv 

In  the  sequel  we  use  freely  the  inequality 
(D  '♦l*2*3,-1  ‘  cl*1l-'*2«0l*3,p  i 

which  is  easily  proved  by  usinq  the  .Sobolev's  inequality  !»lj*  *  cl»l  ,,  Vs  e  Wg*2.  0n* 

has 

(2)  iH(h,v,r,B)  -  H(h,v1,T1,81)l_1  <  l(T-r1)hl_1  +  cvl  (T-r1)vVBI_1  + 

+  I  (m+T1  )(v-v1  )*VBI_1  +  •  (m+T  1  )v1*7  (B-B  ^  )l  _1  + 

+•  l*(v,v)  -  *(v1,v1)l_1  +  M3-ei)(m+T)-1ia2(T,B)(vVT-g)l_1  + 

+  I  (CQtB1  )(m+T  )-1  (m+t  1  )-1ui2(t  ,B  )(vVT-g)l_1  + 

+  l(C0+B1)(m+T1)"1(u)2(T,8)  -  ai2(T1,B1))(vVt-g)I_1  + 

+  •«:0+61)(m+T1)~1ui2(T1,B1)(v-v1)*7Tl_1  + 

+  KC0+61)(m+T1)“1M2(T1,B1)v1*7(T-r1)l_1  . 

By  using  adequately  the  inequality  (1),  one  gets 

(3)  <  clhlplt-r^o  +  ‘’Cy* v*  1#pl®l  1,plx_T  0  + 

+  clm+ljlB^  +  •  (m+r  ^  Jv^V  (8-B  n  )l  + 

+  c(X0+X1)(«vl1p  +  'V^pJIv-v,!,  + 

+  c<2/»>T2('V,1,p,T,1.p  +  I^IP>,S“6 1*0  + 

^(^^^(Wj^dvI^p.T.^p  +  IglpJlT-T,!,,  ♦ 

+  c(q0+l1)(2/m)(lvl1fplTl1>p  v  |,|p)S2(lT-i  4-IS^Ig)  ♦ 

+  ct^+t,  )(2/m)T2lT«1plv-v1l0  + 

+  1  <5„+8  ,)(m+T  )"1w  (t  ,8  )v  *7(t-t  )l 


Me  remark  thee  each  term  on  the  right  hand  side  of  (2)  is  bounded  by  the  term  on  the 

right  hand  side  of  (3)  that  occupies  the  same  relative  position.  In  order  to  estimate  the 

fourth  and  the  last  term  on  the  right  hand  side  of  (3)  we  will  use  the  estimate 

(4)  Iw7*l  <  cdwl.  +  l  |D  w  !  >1*1.  • 

i,j  1  P 

Hare  w  is  a  vector  field  and  a  scalar  field  in  ft.  This  inequality  is  proved  by 
noting  that 

l/n(wV*)fdx!  <  c (  | div  w|p,’M0M2.  ♦  Mjililglfl  ,)  • 

By  using  inequality  (4),  one  easily  verifies  that  the  fourth  and  the  last  term  on  the 
right  hand  side  of  (3)  are  bounded  by  c(mtl«  +  It  ,1  ,  )lv,l  ,  IB-6,1-,  and  by 

1  1  1 1  p  1  1  r  p  1  U 

c(C„+i.  +  18.1.  )(-  +  It.l.  )  (T*  +  S-lx  .1  .  +  8,16.1.  )l  v  l  It  -a  I  , 

0  1  1  1,p  m  2  1  1,p  2  2  1  1,p  2  1  1,p  1  1,p  1  0 

B 

respectively.  By  using  these  estimates  together  with  (3),  one  gets  (3.21). 
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